How quantum bound states bounce and the structure it reveals 
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We investigate how quantum bound states bounce from a hard surface. Our analysis has applica- 
tions to ab initio calculations of nuclear structure and elastic deformation, energy levels of excitons 
in semiconductor quantum dots and wells, and cold atomic few-body systems on optical lattices with 
sharp boundaries. We develop the general theory of elastic reflection for a composite body from 
a hard wall. On the numerical side we present ab initio calculations for the compression of alpha 
particles and universal results for two-body states. On the analytical side we derive a universal 
effective potential that gives the reflection scattering length for shallow two-body states. 



o 

(N 



o 

=5 

> 

oo 

in 

od 
O 

o 



X 



In this letter we consider the elastic scattering of quan- 
tum bound states from a hard surface. Elastic scattering 
from a hard surface provides direct information about the 
structure of the bound state and how it behaves under 
spatial confinement and compression. Spatially-confined 
bound states are interesting quantum systems with im- 
portant practical applications. For example in quantum 
dots and wells, electron- hole pairs are confined to the in- 
terior of nanoscale semiconductor structures. The size 
and shape of the structure strongly influences the density 
of exciton states, thereby allowing control of current tun- 
neling as well as photon absorption and emission. See for 
example Ref . [lHj] ■ In our analysis we present universal 
results for direct band-gap semiconductors, relating ex- 
citon energies to the binding energy at infinite volume, 
the effective masses of the electrons and holes, and the 
geometry of the nanostructure. 

Confined systems can also be produced in cold atomic 
experiments. A quantum well can be constructed us- 
ing an optical lattice with additional lasers focused to 
produce sharp repulsive boundaries. Similar ideas have 
been proposed for quantum billiards systems |5| . Since a 
flat hard-wall boundary has no dimensionful scale, re- 
flection from the wall probes the universal physics of 
few-body systems at large scattering length. For ex- 
ample the dimer-wall reflection scattering length for a 
shallow dimer is directly proportional to the two-particle 
scattering length. In the following we will calculate the 
universal constant of proportionality between these two 
scattering lengths. 

We consider a non-relativistic bound state in d spatial 
dimensions with mass M. The number of constituent 
particles comprising the bound state is arbitrary. The 
bound state scatters elastically against a flat hard-wall 
boundary. Let X be the distance from the wall to the 
center of mass of the bound state. We choose the inertial 
frame where all components of momentum parallel to 
the wall surface are zero and construct a standing wave 
solution with momenta ±p perpendicular to the wall. We 
denote the reflection phase shift as S(p), and define the 
reflection radius as R(jp) = —6{p)/p. This corresponds 
with the distance between the wall and the closest node of 
the asymptotic standing wave, ^ P (X) cx sin [pX + S(p)}. 

The reflection scattering process is analogous to S- 



wave scattering in three dimensions. We can therefore 
use the effective range expansion, 



p cot S(p) =-— + \r R p 2 - Vrp 4 
a R 2 



(1) 



where a r is the scattering length, tr is the effective 
range, and Vr. is the shape parameter. We note 
that aR equals the reflection radius at threshold, clr = 
lirrip^o R{p)- If the bound state structure is completely 
rigid then R(p) — cir for all p. For realistic systems, how- 
ever, compression of the bound state increases with col- 
lision energy. Hence R(p) decreases monotonically with 
p, and the rate of decrease indicates the compressibility 
of the bound state under unilateral stress. The reflec- 
tion radius is directly related to the energy of the bound 
state under spatial confinement. For a d-dimensional 
cube with length L and vanishing Dirichlet boundaries, 
the energy of the lowest confined mode is 



E(L) 



p 2 (L)-d 
2M 
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This formula arises from the sine function profile with 
half-wavelength L — 2R [p(L)} in each dimension. It can 
be used to determine the reflection radius as a function 
of momentum. The O \L~) error estimate contains 
corrections due to double-wall collision effects near the 
wall intersections. 

It is not possible to construct a hard surface for protons 
and neutrons in nuclear physics experiments. Therefore 
it may seem that the current discussion has no direct con- 
nection with nuclear physics. A similar critique could 
be made of Luscher's analysis of periodic boundaries in 
finite cubic volumes [f| 0]- However Luscher's analy- 
sis now provides the theoretical framework for numerous 
calculations in lattice quantum chromodynamics pj-tioj. 
The point is that ab initio numerical calculations play an 
increasingly important role in nuclear theory. In some 
cases the numerical calculations investigate phenomena 
directly observed in experiments. In other cases the 
calculations probe new physics not accessible in experi- 
ments. Similar to temperature and chemical potential, 
different spatial boundaries provide tunable control pa- 
rameters for the numerical calculations. 
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TABLE I: Momenta and reflection radii for an alpha particle 
confined to a cube of length L. 



L 


P(L) 


R[p(L)] 


11.8 fm 


81(9) MeV 


2.1(4) fm 


9.9 fm 


97(10) MeV 


1.6(3) fm 


7.9 fm 


118(10) MeV 


1.3(2) fm 



We now briefly discuss one example from nuclear 
physics. There has been much interest in alp ha-particle 
clusters inside nuclei such as carbon- 12 [lll - [l3| . Very 
recently there have been ah initio lattice calculations of 
the low-lying states of carbon-12 using effective field the- 
ory [3]. In addition to the ground state and excited 
spin-2 state, these calculations find for the first time the 
Hoyle resonance responsible for carbon formation at stel- 
lar temperatures. The lattice data suggests the presence 
of correlated alpha clusters. Furthermore it appears that 
the alpha clusters themselves are compressed in size when 
bound in nuclei. This would be interesting since there 
are no low-energy resonances of the alpha particle and 
little experimental data on the compression of alpha par- 
ticles. 

We have calculated the energy of an alpha particle at 
leading order in chiral effective field theory when con- 
fined to a cubic box with vanishing Dirichlet boundaries. 
We use the same lattice action, algorithms, and codes as 
in Ref. The boundaries are imposed by adding an 

essentially infinite potential energy at the wall bound- 
aries. We note that the ultraviolet divergences are in- 
dependent of long-distance boundary conditions, and so 
no new renormalization counterterms are needed. We 
only summarize the results here, and the numerical de- 
tails of this calculation will be discussed in a forthcoming 
publication [15j . 

At leading order we have calculated the energy for cu- 
bic boxes with lengths L — 11.8 fm, 9.9 fm, and 7.9 fm. 
The results are shown in Table [TJ The error bars in Ta- 
ble [T] are one standard deviation estimates which include 
both Monte Carlo statistical errors and uncertainties due 
to extrapolation at large Euclidean time. The reflection 
radii in Table U can be compared with the root-mean- 
square (RMS) matter radius of the alpha particle. At 
leading order we find the RMS matter radius for pointlikc 
nucleons to be 1.53(4) fm. 

At low momenta the reflection radius of the alpha par- 
ticle is somewhat larger than the RMS matter radius. 
However we see a substantial decrease of the alpha par- 
ticle reflection radius with increasing momentum. This 
data suggests that the alpha particles do compress quite 
readily with pressure. This appears consistent with the 
observation that alpha clusters are compressed in size 
when bound in nuclei. Much more numerical work is 
planned to study alpha particles and other nuclei under 
confinement pressure. 

We now focus on the case of two-body bound states. 
The masses of the two constituent particles will be de- 



noted mi and m2, and fi is the reduced mass, y, = 
{m^ 1 + m^ 1 )^ 1 ■ We assume the two constituent par- 
ticles are distinguishable and so the quantum statistics 
and spin are irrelevant. For the case of electron- hole 
pairs, mi and mi correspond with the effective masses 
of the electrons and holes around zero-momentum min- 
ima. Our analysis applies to the case of Wannier ex- 
citons where the exciton size is significantly larger than 
the Coulomb screening length. 

Let Eb be the energy of the dimer at infinite volume, 
and Kb be the binding momentum defined by the rela- 
tion E B = -k 2 b /(2h). We are assuming that the bound 
dimer has zero orbital angular momentum and let as be 
the S-w&ve scattering length for the two constituent par- 
ticles. The total spin of the bound state is irrelevant 
to our analysis so long as there is negligible partial- wave 
mixing. We consider the shallow-binding limit where 
as is much larger than the range of the interaction. In 
this limit kb = o-b 1 ' ana - the low-energy physics of shal- 
low dimcrs is independent of the short-distance details of 
the interaction. It follows that the dimer- wall reflection 
phase shift is a universal function of the dimensionless 
ratio p/n B . In the following discussion we present the 
form of this universal function. 

We have calculated the dimer-wall reflection phase 
shift for the one- and three-dimensional systems numeri- 
cally using a Hamiltonian lattice formalism. More details 
on this lattice formalism can be found in Ref. [16]. The 
two-particle interaction is implemented as an attractive 
single-site interaction. Two parallel hard-wall bound- 
aries are placed a distance L apart, and we set the mo- 
mentum of the dimer parallel to the wall equal to zero. 
The confinement energy is measured as a function of L us- 
ing sparse-matrix eigenvector methods to determine the 
reflection phase shift. Calculations are repeated sev- 
eral times using successively smaller lattice spacings to 
extrapolate to the continuum limit and determine univer- 
sal results in the shallow-binding limit. For the three- 
dimensional case we also perform an infinite-size extrap- 
olation in the dimensions perpendicular to the wall. A 
large number of separate calculations arc required to per- 
form these extrapolations, and the details of each calcu- 
lation are described in a forthcoming publication (lEj . 

Fig. Q] shows the one-dimensional results for the re- 
flection radius versus dimer kinetic energy, Ek, in the 
shallow-binding limit. The results are presented in di- 
mensionless combinations, kbR(Ek) versus Ek/\E b \. 
We show results for several different mass ratios, mi/mi. 
As expected the reflection radius decreases monotonically 
with energy. One interesting feature is the dependence 
on the mass ratio. For large mi/mi the reflection ra- 
dius is rather large at small energies while significantly 
decreasing with energy. We explain this mass ratio de- 
pendence later in our discussion. 

For the case raijm\ = 1, the one-dimensional problem 
is integrable and can be solved exactly using the Bethe 
Ansatz. From the Bethe Ansatz we get 

pcot(5(p) = ~2kb, (3) 
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m 2 /m i = 1 □ 
m 2 /m, = 2 x 
m 2 /m, = 4 * 
m 2 /m, = 8 + 
Bethe Ansatz m 2 /m 1 = 1 - 
first order Kg3^ O 
second order Kg3^ O 
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FIG. 1: One-dimensional lattice results for the reflection ra- 
dius versus dimer kinetic energy in the shallow-binding limit. 
Also shown are Bethe Ansatz results for 777.2 jm\ — 1 and first- 
and second-order results for the expansion of kbclr described 
later in the text. 

TABLE II: Coefficients of the effective range expansion for 
one-dimensional dimer-wall scattering. 



7772/ 'mi 


KBdR 


KBTR 




1 


0.4999(2) 


0.005(7) 


0.002(3) 


2 


0.6065(2) 


-0.074(2) 


-0.006(2) 


4 


0.8747(2) 


0.115(2) 


0.006(2) 


8 


1.2149(2) 


0.460(2) 


0.008(2) 



As shown in Fig. [TJ the agreement between the lat- 
tice calculations and Bethe Ansatz solution is very good. 
Fig. [U also shows the first- and second-order analytic re- 
sults for the expansion of kbcir described later in our 
discussion. 

The coefficients of the effective range expansion for 
one-dimensional dimer-wall scattering are given in Ta- 
ble [TTJ The error estimates are from least-squares fit- 
ting used in the lattice extrapolation and the effective 
range expansion. As one can see from Eq. ([3]), the Bethe 
Ansatz gives kbcir = 1/2 for mil i m\ = 1, with all other 
higher coefficients equal to zero. This is in full agreement 
with the lattice results. The results are in principle uni- 
versal and can be checked using other theoretical meth- 
ods or perhaps experiments such as cold atomic dimers 
on a one-dimension optical lattice with sharp boundaries. 

Fig. [2] shows the reflection radius versus dimer kinetic 
energy for the three-dimensional system in the shallow- 
binding limit. Again the results are presented in dimen- 
sionless combinations, kbR(Ek) versus Ek/ \Eb\- The 
reflection radius is nearly a factor of two smaller than in 
the one-dimensional case. This is reasonable considering 
the difference between compression of a one-dimensional 
object versus compression of a three-dimensional sphere 



FIG. 2: Three-dimensional lattice results for the reflection ra- 
dius versus dimer kinetic energy in the shallow-binding limit. 
Also shown are first- and second-order results for kbcir de- 
scribed later in the text. 



TABLE III: Coefficients of the effective range expansion for 
three-dimensional dimer-wall scattering. 



1712 /mi 




KBTR 


k%Pr 


1 


0.279(8) 


-0.38(9) 


-0.02(2) 


2 


0.342(8) 


-0.44(6) 


-0.02(2) 


4 


0.506(8) 


-0.27(4) 


-0.02(1) 


8 


0.731(8) 


0.03(2) 


-0.005(5) 



along one direction. As in the one-dimensional case, we 
see the same dependence on the mass ratio. At large 
777,2 /mi the reflection radius is large at small energies 
while becoming substantially smaller with increasing en- 
ergy. Fig.[2]also shows the first- and second-order results 
for the asymptotic expansion of kbcir described later in 
our discussion. 

The coefficients of the effective range expansion for the 
three-dimensional system are given in Table IIIII One 
immediate application of our three-dimensional results 
would be to calculate the reflection phase shift of the 
deuteron from ab initio lattice calculations using chiral 
effective field theory. In addition to verifying the uni- 
versal coefficients in Table IIIII it should be possible to 
measure the spin dependence of the reflection phase shift 
resulting from the D-wave component of the deuteron 
wavefunction. 

We now discuss the analytic calculation of the reflec- 
tion scattering length for a shallow dimer. The calcula- 
tion we present uses the principles of effective field theory 
but covers somewhat unfamiliar territory. Namely, our 
system is inhomogeneous and there is no obvious small 
expansion parameter. 

The method we use is to apply an adiabatic approx- 
imation for the center-of-mass motion and generate an 
asymptotic expansion for the long-distance part of the 
effective potential. We find that this process generates 
an expansion for kbo>r with convergence controlled by 



4 



an expansion parameter of size e _2KBaR . The derivation 
and full details of this expansion up to second order will 
be given in a forthcoming paper [l5| . We summarize the 
main method and results here. 

Let the coordinates of the two constituent particles be 
r*i and r%. The particles have a short-range attractive 
interaction, and we take the Hamiltonian to be 



H 



1 a* 



2m\ 



-VI - 



— V* +CB5^(n-f2), (5) 



2m 2 



where S^ d ' is a regulated d-dimensional delta function. 
The coefficient Cb is tuned to produce a bound state 
with energy Eb at infinite volume. Let r = f± — r 2 be 
the relative separation between the particles. For any 
fixed center-of-mass coordinate the dependence on the 
relative coordinate is given by 



H Icl = -±-V 2 r + C B 6 {d) (r). 



(6) 



To calculate the reflection scattering length it suffices 
to compute dimer-wall scattering in the limit Ek <C 
\Eb\- For this low-energy limit we use an adiabatic ex- 
pansion for the center-of-mass motion. This technique 
is conceptually similar to the adiabatic hyperspherical 
approximation [l7l - [l9| . For each X we keep only the 
ground state for H rc \ satisfying the hard-wall boundary 
condition. Contributions from higher excitations in iJ ro i 
are suppressed by powers of Ek/ \Eb\- Therefore these 
higher excitations contribute to higher-order coefficients 
in the effective range expansion but make no contribution 
to the reflection scattering length. 

For each X let ijjx(r) be the normalized ground state 
wavefunction of H ro i satisfying the hard-wall boundary 
constraint. The low-energy effective Hamiltonian is 



H, 



eff : 



1 d 2 

' 2M dX 2 



+ V(X)+T(X), 



where V(X) is the adiabatic potential, 

V(X) = {^ x \H Iel \^ x ), 
and T(X) is the diagonal adiabatic correction, 



(7) 



(8) 



(9) 



We note that the cross-term involving the expectation 
value of ipx for one derivative with respect to X vanishes 
due to the fixed normalization of ipx ■ 

Let us now define the length scales r+ (X) = 2MX/mi 
and r_ (X) = 2MX/m,2- ?'+(-) PO is twice the distance 
from the wall to ru^} when the other particle touches the 
wall. In other words i/>jc(r) must vanish when the com- 
ponent of r perpendicular to the wall equals r+(X)/2 or 
— r~(X)/2. For large X we can generate an asymptotic 
expansion for V(X) and T(X) in powers of e~ resr +( x ) 
and e -K-sr-(x)_ 

At zeroth order we recover the infinite volume result, 
V(°\X) + TW(X) = E b . At first order the effective 



potential for the one-dimensional case is 
T {1 \X) + V {1) {X) 



k 2 b M 2 
m\m2 



-K B r + (X) 



mi 



3 -KB'r- (X) 



m 2 



(10) 



In three dimensions the result is 



T^(X) + V^(X) 



K B M 

2mim2X 



e -K B r + (X) ^ e -K B r-(X) 



ID 



From these effective potentials it is straightforward to 
compute the first-order result for the reflection scattering 
length. This process can be carried forward to any order. 
The net result is an expansion with an expansion param- 
eter of size e ~ KBr±i - aR ' ) < e - 2K Ba R ^ The larger KbCLr ^ 
the faster the convergence of the expansion. First- and 
second-order results for the one-dimensional system are 
shown in Fig. [TJ and results for the three-dimensional 
system are shown in Fig. [5J In each case the agreement 
with lattice results is consistent with third-order correc- 
tions of size < e ~ 6KBaR . 

We now finally address what happens in the limit 
TO2/TO1 — > 00. Consider the limit m<i 00 with mi 
held fixed. In this limit the effective potential converges 
to a non- vanishing finite-valued function, while the mass 
of the dimer grows with m 2 . One can check this explic- 
itly for the expressions in Eq. (|10l) and Eq. (TTTj) . Given 
the exponential tail of the effective potential, the reflec- 
tion radius near threshold has a logarithmic dependence 
on 77x2/7711. An immediate application of these results 
would be to verify the logarithmic m^/vnx threshold de- 
pendence in various physical systems. This effect should 
be prominent for halo nuclei with a heavy core and one 
satellite nucleon. It should also be seen in quantum dots 
and wells for semiconductors with a large ratio between 
hole and electron effective masses. It could also be repro- 
duced with heterogeneous cold atomic dimers consisting 
of one heavy and one light alkali atom. 

In this letter we have discussed a number of phenom- 
ena relating to the elastic scattering of quantum bound 
states from a hard surface. We have emphasized pro- 
cesses of universal character which can be seen in several 
different systems. The applications appear numerous, 
ranging from experimental predictions for quantum dots 
and wells to numerical calculations of nuclear structure. 
In our analysis we have introduced some theoretical tech- 
niques which may be useful for describing the effective 
field theory of other inhomogeneous systems. 
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